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I  INTRODUCTION  AND  SUMMARY 


A.  Motivation  and  Objectives 

The  determination  of  dynamic  stresses  and  loads  is  fundamental  to 

much  of  DNA  field  testing.  The  need  for  these  measurements  has  led  to 

the  development  and  use  of  many  different  types  of  stress,  particle 
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velocity,  acceleration,  and  displacement  gages.  ’  Despite  the  large 
variety  and  quantity  of  existing  dynamic  measurements,  field  measurements 
of  shear  stress  (or  loads)  are  lacking.  The  inability  to  make  dynamic 
shear  measurements  is  an  important  shortcoming,  because  shear  measure¬ 
ments  are  needed  for  determining  strength  properties  of  both  soils  and 
structures  in  underground  tests. 

The  need  for  shear  measurements  is  well  recognized  by  most  workers 
involved  in  dynamic  measurements.  However,  this  development  has  been 
lacking  due  to  the  complexity  of  the  problem:  The  desired  shear  stress 
gage  must  be  usable  under  complex  loading  conditions,  and  suitable 
methods  are  needed  to  calibrate  the  gage  to  well-defined  shear  stresses. 
Most  of  the  previous  field  gage  techniques  have  been  extended  from 
laboratory  concepts  and  measurements,  but  laboratory  studies  are  also 
lacking  in  dynamic  shear  measurements. 

The  objective  of  our  work  was  thus  to  examine  the  feasibility  of 

developing  a  piezoelectric  shear  stress  gage  for  use  in  DNA  field  tests. 

A  combined  analytic  and  laboratory  experimental  effort  was  undertaken 

to  meet  this  objective.  The  bases  for  this  work  are  recent  developments 

3  4 

at  SRI  relating  to  the  study  of  dynamic  shear  properties  of  solids 

and  the  use  of  piezoelectric  gages  in  studying  dynamic  compressive 

5-9 

stresses. 


B.  Background  and  Approach 


Piezoelectricity  and  the  use  of  piezoelectric  transducers  and 

devices  in  acoustic  applications  is  a  major  field  of  study. ^  ^ 

However,  the  use  of  piezoelectric  transducers  for  measuring  large 

stress  amplitudes  under  dynamic  loading  is  a  more  recent  and  specialized 
13 

topic.  Studies  during  the  past  decade  have  led  to  the  development  of 
stress  transducers  for  studying  dynamic  compressive  stresses  up  to  tens 
of  kilobars  in  materials  subjected  to  impact,  explosive,  and  radiation 


loads. 
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Studies  by  Graham  and  co-workers  have  clarified  piezo¬ 


electric  response  at  high  stresses  and  established  bounds  on  the  use  of 
stress  gages. ^  ^  The  formulation  of  nonlinear  piezoelectric  con¬ 
stitutive  relations  has  also  received  increased  attention. In 
many  laboratory  and  field  measurements,  alpha  quartz  is  used  as  the 
gage  material,  but  recent  studies  have  also  been  conducted  using 

TO 

lithium  niobate  (LiNbO-^).  For  low  stresses  (below  10  kbar),  the 
larger  electrical  output  of  LiNbO-  is  advantageous. 


Piezoelectric  gages  are  used  in  two  modes.  The  current  or  short- 
circuit  mode  measures  fast  ri  .e,  short-duration  stress  pulses, 

and  the  useful  recording  time  o  the  gage  is  the  wave  transit  time 
through  the  gage.  Laboratory  shock  wave  experiments  with  zero  lateral 
strains  commonly  use  the  current  mode.  In  the  charge  or  open-circuit 
mode,  the  gage  is  used  to  record  slower  rise  time,  long-duration 
(milliseconds)  pulses.  This  second  mode  is  more  commonly  used  under 
field  conditions  for  recording  stress  pulses  with  wavelengths  much 
greater  than  the  gage  thickness.  That  is,  the  gage  in  the  charge  mode 
acts  like  a  static  transducer  in  equilibrium  with  the  surrounding 
material  stresses.  No  fundamental  differences  exist  between  these 
modes  \\?ith  regard  to  the  piezoelectric  response,  and  the  laboratory 
results  are  applicable  to  field  usage.  The  differences  in  these 
modes  are  operational  and  reflect  the  mechanical  boundary  conditions 
most  suited  for  using  the  gages  in  different  applications.  Further 
discussion  of  these  two  operational  modes  is  given  in  Reference  6  and 
Appendix  A. 
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Despite  the  many  laboratory  and  field  developments  in  the  use  of 
the  above-cited  modes,  measurements  to  date  have  been  performed 
exclusively  for  compressive  stresses-  In  the  following  paragraphs  we 
discuss  our  approach  for  developing  a  piezoelectric  shear  stress  gage 
for  use  in  dynamic  loading  conditions  in  the  presence  of  complex 
stresses. 

The  two  general  requirements  in  the  uevelopment  of  a  shear  gage 
are:  (a)  knowledge  of  a  material  phenomenon  (e.g  ,  appropriate 

piezoelectric  response)  relating  shear  stress  to  a  measurable  quantity 
and  (b)  the  ability  to  calibrate  the  gage  for  known  shear  stresses. 

The  first  requirement  is  easily  satisfied,  in  principle,  because 
of  the  existence  of  many  shear  transducers  in  the  field  of  ultrasonics. 
Because  the  compressive  stress  gage  is  an  extension  of  ultrasonic 
concepts,  these  concepts  can  be  explored  for  the  development  of  a  shear 
stress  gage.  There  are,  however,  complicating  factors  that  do  not 
allow  a  simple  extension  of  ultrasonic  concepts.  In  ultrasonics, 
pure  shear  waves  are  commonly  created  by  using  the  converse  piezo¬ 
electric  effect.  Thus,  the  measuring  transducer  is  subjected  to  a  pure 
shear  motion.  In  most  dynamic  loading  situations,  a  complex  stress 
state  exists — that  is,  superposed  compression  and  shear  states. 
Furthermore,  because  of  the  large  stresses,  the  material  in  which  the 
gage  is  placed  undergoes  inelastic  deformations  and  the  relative 
magnitudes  of  the  compression  and  shear  stresses  are  expecteo  to  vary 
over  the  time  range  of  interest.  For  the  gage  measurements  to  be 
useful,  the  electrical  signal  from  the  gage  should  be  uniquely  related 
only  to  the  shear  stress  of  interest.  Because  of  the  tensorial  nature 
of  piezoelectricity,  this  requirement  is  not  met  by  most  of  the 
ultrasonic  shear  transducers. 

To  satisfy  the  above  requirement  in  the  presence  of  complex 
stresses  (e.g.,  combined  compression  and  shear),  we  developed  an 
analytic  approach  that  consists  of  simultaneously  examining  wave 
propagation  and  piezoelectric  relations  for  materials  of  interest. 

The  analytic  approach  chosen  is  quite  general  and  provides  criteria 


for  determining  the  needed  gage  designs.  We  examine  all  possible 
crystal  orientations  to  obtain  the  optimal  directions  for  piezoelectric 
response  and  to  determine  pure  mode  directions  ('specific  directions') 
for  stress  wave  propagation.  Once  the  optimal  gage  designs  have  been 
analytically  evaluated,  we  can  conduct  experiments  to  verify  and 
calibrate  these  gage  designs. 

The  second  requirement,  calibration  of  the  gages,  is  met  by 

experimentally  studying  the  gage  response  to  one-dimensional  compression 

* 

and  shear  waves  using  the  recently  developed  IMPS  experimental 

4+ 

facility.  ‘  In  this  method,  the  specimens  are  subjected  to  varying 
but  controlled  amounts  of  compression  and  shear  stresses.  By  subjecting 
the  gages  to  pure  compression,  we  can  ensure  that  there  is  no  electrical 
output  from  the  compressive  stress.  Increased  amounts  of  shear  stresses 
can  then  be  superposed  to  calibrate  the  response  to  shear.  Because  this 
is  a  feasibility  study,  the  present  scope  of  the  work  is  intended  to 
experimentally  verify  the  theoretical  concepts.  A  detailed  calibration 
would  be  performed  after  the  feasibility  of  the  shear  gage  has  been 
established. 

C.  Summary 

The  objective  of  the  work  reported  here  was  to  examine  the  feasi¬ 
bility  of  developing  a  piezoelectric  shear  stress  gage  sensitive  only 
to  shear  loading.  Using  a  combined  analytic  and  experimental  approach, 
we  successfully  demonstrated  the  feasibility  of  developing  such  a  gage. 

A  simplified  but  general  analysis  of  the  mechanical  and  piezo¬ 
electric  response  was  carried  out  to  provide  criteria  for  a  suitable 
shear  gage.  Mechanical  wave  propagation  analysis  in  anisotropic 
electic  media  showed  that,  in  general,  three  waves  are  propagated: 
one  quasilongitudinal  and  two  quasitransverse  waves.  Only  for  'specific 

directions'  are  the  propagated  waves  purely  longitudinal  and/or  purely 
■k 
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transverse.  To  avoid  mechanical  coupling  of  strains,  the  gage  thickness 
direction  should  be  along  a  specific  direction.  The  piezoelectric 
analysis  showed  that  electrical  polarization  for  the  desired  gage 
should  be  one-dimensional  and  along  the  gage  thickness  direction. 
Furthermore,  this  electrical  polarization  should  be  caused  only  by 
shear  loads  and  not  compressive  loads. 

The  adequacy  of  a  particular  crystal  type  for  meeting  the  above 
criteria  can  be  easily  and  efficiently  assessed  by  numerical  calcu¬ 
lations  of  the  specific  directions  and  the  'piezoelectric  matrix'  for 
all  possible  orientations  about  the  three  crystallographic  axes. 

When  these  calculations  were  carried  out  for  the  different  orientations 
of  alpha  quartz  and  lithium  niobate  (LiNbO^),  only  the  163°  Y-cut 
LiNbO^  appeared  suitable  for  use  as  a  shear  gage.  For  this  orientation, 
the  shear-to-compression  sensitivity  ratio  for  the  polarization  is  more 
than  two  orders  of  magnitude.  Furthermore,  this  orientation  deviates 
only  2°  from  a  specific  direction  and  therefore  has  minimal  mechanical 
coupling. 

Impact  experiments  were  conducted  under  combined  compression  and 
shear  to  verify  the  use  of  163°  Y-cut  LiNbO^  as  a  shear  stress  gage. 

The  results  of  these  experimen  -■  show  that  the  gage,  as  desired,  had 
negligible  sensitivity  to  compressive  stress  and  a  very  large  sen¬ 
sitivity  to  shear  stress  (Section  III). 

Further  development  of  gages  for  routine  use  in  laboratory  and 
field  measurements  requires  caliDration  experiments  that  can  be  per¬ 
formed  uring  the  same  impact  facility.  In  addition,  field  usage 
requires  development  of  proper  packaging  techniques.  A  brief  dis¬ 
cussion  of  these  calibration  and  packaging  requirements  is  also 
presented  in  this  report  (Section  IV) . 


II  THEORETCIAL  CONSIDERATIONS 


This  section  describes  the  bases  for  selecting  the  gage  designs 
used  in  our  impact  experiments.  We  present  the  simplifying  assumptions 
made  in  our  analysis,  followed  by  an  examination  of  elastic  wave 
propagation  and  piezoelectric  relations  for  the  materials  of  interest. 
Then  we  describe  the  specific  gage  designs  selected  on  the  basis  of 
these  analytic  concepts.  The  analysis  presented  here  is  not  concep¬ 
tually  new,  but  rather  is  a  systhesis  of  existing  studies  ’ *  ’ ^  for 
our  objectives. 

We  are  primarily  concerned  with  two  specific,  crystals:  alpha- 
quartz  and  lithium  niobate  (LiNbO^).  These  crystals  have  been  used 
in  previous  dynamic  measurements  because  of  their  high  mechanical 
strength  and  easy  availability. 

A.  Simplifying  Assumptions 

Finite  strains  and  large  electric  fields  are  created  in  the  piezo¬ 
electric  gages  subjected  to  large  amplitude  stress  waves.  Thus,  a 
detailed  discussion  of  the  piezoelectric  constitutive  relations  should 

incorporate  finite  strain  effects,  electromechanical  coupling,  and 

15 

changes  in  d! electric  permittivity.  Analysis  of  the  experimental 
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results  should  also  consider  these  nonlinear  effects,  ’  but  such 
an  effort  is  beyond  the  scope  of  a  feasibility  study. 

In  the  present  work,  we  make  use  of  two  simplifying  assumptions: 
"linearity"  and  "no  coupling. The  first  assumption  gives  a  linear 
stress-strain  relation  and  a  linear  piezoelectric  relation.  The  second 
assumption  permits  us  to  neglect  the  effect  of  electric  fields  on 
mechanical  variables  such  as  strains  and  wave  velocities.  Therefore, 
we  can  tudy  the  piezoelectric  response  in  two  separate  parts:  elastic 
wave  propagation  (mechanical  loading)  and  piezoelectric  relations. 

Each  of  these  two  parts  suggest  desirable  features  for  the  gage 
design. 


B.  Elastic  Have  Propagation 

Wave  propagation  in  a  linear  elastic  anisotropic  medium  has  been 

19-21 

the  subject  of  considerable  work  in  ultrasonics  and  acoustics. 

Johnson  has  applied  these  concepts  to  planar  impact  problems  and 

has  obtained  analytic  solutions  that  show  multiple  wave  structure 

22 

because  of  crystalline  anisotropy.  Below,  we  summarize  the  main 

* 

features  of  elastic  wave  propagation  in  anisotropic  crystals. 

22 

We  consider  two  coordinate  systems:  x^  corresponds  to  the 

crystallographic  system  and  x'  corresponds  to  the  wave  propagation 

1  20 

system.  Wave  propagation  is  described  by 
32u:  32u" 

P  ~  2  =  Cijk£  3x~:  3x'"  (1) 

3t  1  l 


where  p  =  density 

u'  =  material  displacement  along  the  x7-direction 
t  =  time 

=  second-order  elastic  constants  in  the  primed  system. 

The  transformation  matrix  between  the  unprimed  and  primed  system  is 
given  by 

~>v  ->■ 

a  =  i  •  i  (2 

mn  m  n 

—y  ^  ->■  A 

where  i  and  i  are  unit  vectors  along  the  x  and  x  directions, 
m  n  °  m  n 

-V 

Plane  wave  solutions  for  wave  propagation  along  a  direction  d 
for  Equation  (1)  are  given  by 


uT  *=  u:  f(t  -  b'x'/c') 

i  x  n  n 


where  IK  =  displacement  amplitude  along  xT 
c'  =  wave  velocity. 

Substitution  of  (3)  in  (1)  gives  the  following  equations 


The  Einstein  summation  convention  is  adopted  throughout  this  report. 
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Aik  -  sik  pc' 


where  A'k  =  b'h'  C1Jtl 
For  nontrivial  solutions  of  Eq.  (4),  we  have  the  conditon 


iA;k-5ik  pc  i  -  °  «> 

Equation  (6)  has  three  possible  eigenvalues,  each  of  which  corresponds 

to  an  eigenvector  U^.  Each  of  the  eigenvalues  and  eigenvectors 

represent  the  wave  velocity  and  displacement  amplitude  associated  with 

one  of  the  three  waves.  These  displacements  are  always  mutually 

orthogonal  but  can  have  an  arbitrary  orientation  with  respect  to  the 

wave  propagation  direction  d  .  Only  for  specific  directions  (or 

isotropic  materials)  are  the  displacement  directions  either  parallel 

(longitudinal  waves)  or  perpendicular  (shear  waves)  to  the  direction 
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of  wave  propagation.  In  general,  there  is  one  quasilongitudinal 
wave  and  two  quasitransverse  waves. 

In  our  experiments  we  are  subjecting  the  gage  to  externally  applied 
compressive  and  shear  stresses.  If  the  gage  output  is  to  be  related 
to  a  particular  stress  component  (shear  stress  in  the  present  case), 
then  no  coupling  of  the  stresses  should  occur  within  the  gage.  For 
the  linear  elastic  behavior  considered  here,  this  coupling  can  arise 
only  because  of  the  crystal  anisotropy.  Therefore,  the  crystals  used 

"k 

for  gage  development  should  be  oriented  along  the  specific  directions, 
for  the  compressive  and  shear  wave  of  interest  in  our  experiments. 

The  specific  directions  for  a  particular  crystal  depend  on  the 
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symmetry  of  that  crystal  class.  Borgnis  and  Brugger  have  performed 
extensive  calculations  of  specific  directions  for  various  crystal 
classes.  Because  in  practical  applications  small  deviations  from 
specific  directions  can  be  tolerated,  we  have  chosen  to  examine  wave 
k 

Specific  directions  can  vary  depending  on  the  wave  motion  of  interest 
(the  longitudinal  or  either  of  the  transverse  waves). 


Both  LiNbO,,  and  quartz  are  trigonal  crystals. 


propagation  for  quartz  and  LiNbO^  more  generally.  This  is  best  done 
by  solving  Eq.  (6)  to  give  the  eigenvalues  and  eigenvectors  for  each 
one  degree  of  rotation  in  the  XY,  YZ,  and  ZX  planes  about  the  Z,  X, 
and  Y  axes,  respectively.  We  have  numerically  carried  out  this  cal¬ 
culation  for  the  YZ  plane  (rotation  about  X  axis)  and  the  results 
are  discussed  in  Section  II. D.  The  calculations  for  XY  and  ZX  planes 
were  not  performed  because  of  piezoelectric  considerations  discussed 
in  the  next  two  subsections. 


C.  Piezoelectric  Relations 

The  direct  piezoelectric  effect  for  the  linear  approximations 
made  in  our  work  may  be  expressed  as"*^ 


?' 

1 


el.,  £' 

ijk  jk 


(7) 


where  P'  =  electrical  polarization  along  xT 

eijK  =  piezoelectric  stress  coefficient  tensor  in  the  primed 
system 

=  strain  coefficients  in  the  primed  system. 

For  convenience  in  our  impact  experiments,  we  have  chosen  to  express 
the  piezoelectric  effect  in  terms  of  mechanical  strains  instead  of 
stresses.  On  the  other  hand,  for  the  field  measurements,  it  is 
convenient  to  express  the  polarization  in  terms  of  the  stresses. 


PT  =  dT  o' 
l  imn  mn 


(8) 


where  d'  =  piezoelectric  strain  coefficent  tensor  in  the  primed 
imn  1  1 

system 

=  stresses  in  the  primed  system. 

Using  Hooke's  Law,  we  can  express  the  relationship  between  the  d's 
and  e's  as 


dT  C' 
xpq  pqjk 


(9) 
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We  are  ignoring  pyroelectric  effects  in  our  discussion. 


Equations  (7)  and  (8)  are  equivalent  representations,  and  the  use  of 
either  form  is  a  matter  of  convenience  based  on  the  mechanical  boundary 
conditions. 

Equation  (7)  can  be  expressed  in  matrix  notation  as"^ 


where  we  have  replaced  the  strain  matrix  z'  by  the  column  vector  z' , 

jk  '  n’ 

and  the  off  diagonal  elements  of  are  given  by  n>  3.  For  the 
desired  gage,  we  have  the  following  requirements:  (1)  PT  should 
have  only  one  component,  (2)  this  component  should  be  produced  only 
by  a  shear  strain  z'  (where  n>3),  and  (3)  the  component  should  be 
parallel  to  the  gage  thickness  direction. 


We  now  examine  the  piezo] ectric  response  for  the  loading  situation 
present  in  our  impact  experiments.  Let  X'  be  the  direction  through  the 
gage  thickness  and  z'~  and  €.'  be  the  non-zero  strains.  Calculation  of 
polarization  solely  caused  by  z'^  requires  that  the  piezoelectric 
matrix  e7^  in  Equation  (10)  possess  a  particular  form.  To  determine 
the  adequacy  of  alpha  quartz  (crystal  class  32)  and  LiNbO^  (crystal 
class  3m)  for  our  purposes,  we  first  write  down  the  piezoelectric 
tensor  in  the  crystallographic  system  (X^  system) ^ 


Quartz: 


0 

0 


e14 

0 


0  0  0 


0 

-e14 

0 


(ID 


Because  the  gage  behaves  elastically,  either  stresses  or  strains 
can  be  used  in  the  discussion. 


14 


/ 0 

0 

0 

0 

e15 

"2e22\ 

("e22 

e22 

0 

e15 

0 

0  ' 

1  (12) 

\  e31 

e31 

e33 

0 

0 

°  ) 

LiNbO^ : 


Of  the  three  directions  in  the  crystallographic  system  (X^  system), 
only  Y-cut  quartz  and  X-cut  and  Y-cut  LiNbO^  will  give  non-zero  polari¬ 
zation  for  applied  shear  strains.  Therefore,  we  have  the  following 
situation  for  combined  compression  and  shear  loading  (present  in  our 
experiments)  for  these  orientations: 


Y-cut  quartz: 

(e2  7*  0,  e6  ^  0) 

X-cut  LiNbO^: 

(ex  t  0,  e6  +  0) 

Y-cut  LiNbO^: 

(e2  $  0,  /  0) 


P1  ^  0,  p2  j  0 


P,  *  0,  P,  t  0,  P_  *  0 


PjL  4  0,  p2  4  0,  p3  /  0 


In  the  presence  of  pure  shear  strains  (as  is  often  the  case  in  ultra¬ 
sonic  applications),  the  above  orientations  are  satisfactory  because 
only  one  polarization  term  will  be  non-zero.  However,  in  the  presence 
of  compressive  loading,  these  orientations  give  rise  to  unwanted 
polarizations.  Thus,  for  complex  loading  situations  we  must  include 

the  effect  of  additional  polarization  fields  and  not  just  the  polari- 
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zation  field  for  shear.  ’  ’  “  A  possible  method  to  minimize  the 

effect  of  these  unwanted  polarization  fields  is  discussed  in  the 
next  subsection. 

A  fundamental  method  of  eliminating  the  undesired  polarization 
fields  is  to  find  a  suitable  crystal  orientation  such  that  the  piezo¬ 
electric  tensor  has  the  desired  form — that  is,  where  the  polarization 
field  is  one-dimensional  and  is  caused  solely  by  shear  loads.  The 
transformation  for  the  piezoelectric  tensor  is  expressed  as 


Gijk  ai£ajmakne£mm  ’ 
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An  efficient  way  to  determine  the  crystal  directions  of  interest  is 
to  numerically  transform  the  tensor  for  each  one-degree  rotation 

in  the  XY-plane  (rotation  about  the  Z-axis),  the  YZ-plane  (rotation 
a- Jut  the  X-axis),  and  the  ZX-plane  (rotation  about  the  Y-axis).  By 
expressing  the  transformed  piezoelectric  constants  in  a  matrix 
notation  similar  to  that  of  Equation  (11)  or  (12),  we  can  quickly 
scan  through  the  results  and  determine  the  desired  orientation.  The 
numerical  method  and  results  are  presented  in  Appendix  B.  We  have 
used  this  method  to  examine  the  piezoelectric  tensor  for  alpha  quartz 
and  LiNbO^,  and  the  results  of  these  calculations  are  discussed  below. 

D.  Optimal  Gage  Design 

The  development  of  an  optimal  gage  design  has  three  requirements 
based  on  theoretical  considerations:  (a)  the  electric  polarization 
must  be  due  to  shear  stress  (or  strain)  in  the  gage,  (b)  the  polari¬ 
zation  must  be  one-dimensional  and  along  the  gage  thickness  direction 
and,  (c)  the  gage  thickness  direction  must  be  a  specific  direction. 

We  have  examined  alpha  quartz  and  LiNbO^  in  terms  of  these  requirements. 
Appendix  B  presents  the  results  of  the  transformed  piezoelectric 
matrices  for  rotation  about  the  X-axis  for  LiNbO^.  The  transformed 
matrices  for  quartz  and  the  other  two  axes  in  LiNbO^,  and  calculations 
of  specific  direction,  have  also  been  carried  out  but  are  not  presented 
here. 

In  alpha  quartz  it  is  not  possible  to  eliminate  the  unwanted 
polarizations  by  using  a  rotated  crystal  cut.  The  possibility  of 
mitigating  the  effects  of  such  polarizations  by  suitable  circuit 
design  is  considered  later  in  our  discussion  on  LiNbO^. 

The  results  of  the  calculation  for  LiNbCL  show  that  the  163° 

* 

Y-cut  is  the  optimal  direction  from  polarization  considerations. 

The  numerical  values  of  the  piezoelectric  matrix  (eC.)  for  this 

2 

orientation  expressed  in  the  units  of  Coulomb/mt  are: 


163°  Y-cut  implies  a  163°  rotation  about  the  X-axis.  See  Institute 
of  Radio  Engineers  (IRE)  convention  (Reference  23). 
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.39  -8  x  10 


-1.72  -3.76 


-2.27  -  .69 


The  transformed  values  show  that  for  e'  /  0  and  #  0  (the  condition 
of  interest  in  our  impact  experiments),  the  shear  contribution  to 
is  more  than  two  orders  of  magnitude  larger  than  the  compressive 
contribution.  Although  is  zero,  for  this  orientation,  the  value 
of  P^,  although  small,  is  not  completely  negigible.  The  need  to 
ascertain  the  influence  of  P'  is  discussed  in  Section  IV. 

The  results  presented  in  Appendix  B  show  the  marked  sensitivity 
of  the  piezoelectric  tensor  to  small  changes  in  the  angles.  Both  the 
162°  and  164°  Y-cut  directions  have  lower  shear-to-compression 
sensitivity  ratios  in  comparison  to  the  163°  angle.  These  results 
emphasize  the  need  for  accurate  knowledge  of  the  piezoelectric  constant 
(e^k  in  the  crystallographic  system)  and  for  precise  crystal  orienta¬ 
tion  in  gage  fabrication. 

The  163°  Y-cut  direction  is  also  a  desirable  orientation  on  the 
basis  of  the  mechanical  considerations.  We  have  computed  the  eigen¬ 
values  and  eigenvectors  of  Equation  (6)  for  all  directions  in  the 
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Y-Z  plane,  using  constant  electric  displacement  elastic  constants. 

These  calculations  show  that  the  165°  Y-cut  orientation  is  a  specific 
direction.  Also,  the  longitudinal  particle  motion  direction  for 
wave  propagation  along  the  163°  Y-cut  orientation  is  along  161.7  , 
and  hence  the  mechanical  coupling  effects  due  to  crystal  anisotropy 
are  quite  minimal.  The  complete  se.  of  eigenvectors  for  the  163° 
orientation  is:  (1,  0,  0),  (0,  -0.949,  0.3145),  and  (0,  0.3145, 

0.949). 

The  use  of  163°  Y-cut  LiNbO,  for  acoustic  applications  has  been 

2  A 

proposed  by  Warner,  et  al.  However,  the  criteria  derived  in  our  work 
& 

are  more  general  and  can  be  used  to  evaluate  any  crystal  for  use  as 
a  shear  stress  gage  under  dynamic  loading. 

Vor  example,  our  analysis  shows  that  alpha  quartz  is  not  suitable 
as  a  shear  gage. 


These  analytic  considerations  show  that  even  the  optimal  gage 
direction  for  shear  response  is  not  ideal,  as  is  frequently  the  case 
for  compressive  stress  gages.  Further  discussion  of  the  departures 
from  ideal  considerations  is  given  in  Section  V. 

We  also  conducted  experiments  on  X-cut  LiNbO^,  because,  some 

X-cut  crystals  were  immediately  available  and  because  we  wanted  to 
examine  a  possible  method  of  minimizing  the  effect  of  the  unwanted 

polarizations.  We  wanted  to  find  out  whether  the  'shorted'  gage 
configuration  (the  gage  sides  are  also  plated)  would  aid  in  minimizing 
the  effects  of  the  unwanted  polarization. 


Ill  EXPERIMENTAL  METHOD  AND  RESULTS 


The  gage  designs  determined  from  theoretical  considerations  were 
experimentally  checked  by  performing  impact  experiments  under  combined 
compression  and  shear  loading.  This  section  describes  the  gage 
fabrication,  experimental  techniques  used,  and  the  results  of  these 
impact  experiments. 

A.  Experimental  Design  and  Method 

All  of  the  LiNbO^  gages  used  in  our  work  were  fabricated  by 
Specialty  Engineering  Associates  according  to  our  specifications. 

The  LiNbO.,  crystals  were  cut  to  the  desired  orientation  (within  one 
degree).  The  flatness  and  parallelism  requirements  were  typical  of 
piezoelectric  gages  used  in  impact  experiments.  All  the  gages  were 
gold-plated  on  all  sides  (i.e.,  wrap  around  or  shorted  configuration). 
Guard  rings  were  cut  on  the  back  side  and  a  50-ohm  resistor  was  shunted 
across  the  guard  ring.  The  gages  were  built  into  an  aluminum  housing 
with  epoxy  potting  at  the  rear  of  the  gage.  These  gages  were  designed 
to  record  useful  information  for  one-wave  transit  time  through  the 

gage.  The  guard-ring  and  gage  dimensions  were  similar  to  those  used 
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in  studies  on  compressive  stress  gages. 

All  the  experiments  performed  in  this  work  employed  the  back 
surface  gage  configuration  shown  in  Figure  1.  A  flyer  plate  inclined 
at  an  angle  90°  -0  to  the  axis  of  the  projectile  motion  impacts  a 
parallel  specimen  plate  backed  by  the  piezoelectric  gage.  The  shear 
axis  on  the  gage  plane  has  to  be  accurately  aligned  to  the  direction 
of  shear  stress  imparted  by  the  projectile  impact.  Upon  impact, 
two  waves  are  propagated  into  the  flyer  and  specimen:  a  compressive 
and  a  shear  wave.  The  compressive  wave,  which  travels  at  a  faster 
velocity,  impacts  the  gage,  resulting  in  a  transmitted  wave  into  the 

A  , 

Future  gages  will  be  fabricated  with  tolerances  to  within  ±  1/2°. 


gage  and  a  reflected  wave  in  the  specimen.  The  slower  traveling  shear 

* 


wave  interacts  with  the  reflected  compressive  wave  and  arrives  at 
the  specimen-gage  interface.  For  a  gage  responding  only  to  shear, 
we  expect  to  see  no  signal  from  the  compressive  wave  and  a  large 
signal  after  the  arrival  of  the  shear  wave. 


The  use  of  a  back  surface  gage  for  determining  gage  calibrations 
is  complicated  because  it  is  necessary  to  have  accurate  knowledge  of 
the  shear  response  of  the  specimen,  the  wave  interactions  in  the 
specimen,  and  the  shear  strength  of  specimen-gage  interface.  However, 
the  accurate  determination  of  gage  calibration  was  not  the  objective 
of  this  study.  Instead,  we  intended  only  to  examine  the  experimental 
feasibility  of  the  theoretically  determined  gage  designs.  The  back 
surface  gage  configuration,  because  of  the  temporal  separation  between 
the  compressive  and  shear  waves,  permits  an  efficient  assessment  in 
a  single  experiment  of  the  desired  gage  features:  no  signal  from  the 
compressive  wave  and  a  large  signal  from  the  shear  wave.  In  future 


experiments  for  determining  quantitative  calibrations,  we  will  use 
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the  impact  surface  gage  configuration. 


The  impact  experiments  were  performed  using  the  SRI  gas-gun 
facility  for  studying  one-dimensional  compression  and  shear  wave 
propagation  in  solids.  Details  of  the  gun-barrel,  target  holder,  and 
related  instrumentation  for  measuring  projectile  velocity  and  impact 
misalignment  are  given  in  Reference  4,  Seven  experiments  were  performed: 
two  experiments  with  X-cut  LiNbO-j  gages  and  five  with  163°  Y-cut 
LiNbO^  gages.  The  LiNbO^  gages  were  bonded  to  the  specimen  plates 
(Polymethyl-methacrylate  discs)  using  epoxy  resin  (Epon  815  and 
hardener').  The  specimen-gage  assemblies  were  potted  in  aluminum 
target  rings  as  shown  in  Figure  2.  In  some  experiments,  front  surface 
particle  velocity  gages  were  vapour-deposited  on  the  PMMA  specimen 


In  our  present  work,  we  assume  linearity  and  have  ignored  this 
interaction. 


1  Shell  Co. 
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discs  to  record  the  particle  velocities  imparted  to  the  specimens. 

Tilt  pins  were  vapour-deposited,  and  trigger  foils  were  bonded  to  all 
the  specimens  as  shown  in  Figure  2.  The  flyer  plates  for  the  two  X-cut 
LiNbO^  and  the  three  163°  Y-cut  I.iNbO^  experiments  were  PMMA  discs. 

The  remaining  two  experiments  used  flyer  plates  made  from  syntactic 
foam  ECCO  float  EF-38.  The  rationale  for  using  the  foam  flyer 
plates  is  discussed  in  the  next  section. 

B.  Experimental  Results 

Table  I  summarizes  the  experimental  details  from  the  seven  impact 
experiments  conducted.  The  first  nine  columns  are  self-explanatory. 

The  longitudinal  compressive  stresses  cited  in  the  last  column  are 
the  stresses  in  the  LiNbO^  gage  based  on  a  linear  elastic  calculation: 


2Z 


JL 


zL  + zL 

g  P 


(15) 


where 


O  =  Longitudinal  stress  in  the  gage 


Longitudinal  mechanical  impedance  of  the 
gage  (function  of  crystal  orientation) 


Longitudinal  mechanical  impedance  of 
PMMA  specimen 


a  =  Longitudinal  stress  caused  by  the  impact 
^  in  PMMA 

The  longitudinal  particle  velocities  in  the  PMMA  are  either  measured 

or  known  from  the  symmetric  impact  conditions.  These  velocities  in 

turn  permit  a  determination  of  Z^  and  0^  using  the  data  of  Barker  and 
27  P  P  L 

Hollenbach.  The  mechanical  impedance  of  the  LiNbO^  (Z^)  is  known 

from  the  wave  velocity  calculations  outlined  in  the  last  section. 


Trademark,  Emerson  and  Cuming  Inc.,  Floatation  Products  Division, 
Canton,  Mass. 
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EXPERIMENTAL  PARAMETERS  FOR  LiNb0o  IMPACT  EXPERIMENTS 


The  first  series  of  experiments  were  conducted  on  X-cut  LiNbO^. 

The  gage  dimensions  used  were  the  sizes  that  were  readily  available. 

The  first  shot  was  a  uniaxial  strain  shot  to  determine  the  effect  of 
the  shorting  configuration  in  reducing  the  polarization  charges  per¬ 
pendicular  to  the  direction  of  interest.  Figure  3(a)  shows  the  voltage¬ 
time  profile  from  the  shot  normalized  with  respect  to  the  transit 
time  of  the  signal  through  the  gage.  The  observed  signal  is  very  large 
and  clearly  demonstrates  the  inability  of  the  shorting  configuration 
to  eliminate  the  unwanted  polarization  charges. 

A  second  X-cut  LiNbO^  shot  was  performed  with  the  same  compressive 
stress  as  the  first  shot  but  with  a  superposed  shear  stress.  The 
voltage-time  profile  from  this  shot  is  shown  in  Figure  3(b).  An 
algebraic  difference  between  these  two  shots  is  shown  in  Figure  3(c). 
This  difference  can  be  attributed  to  shear  loading.  Although  the  signal 
in  Figure  3(c)  is  quite  large  (approximately  10  volts),  no  quantitative 
interpretation  of  these  results  is  attempted  here.  The  presence  of 
?2  and  P'  polarization  terms  do  not  permit  the  use  of  assumptions  needed 
to  derive  the  current-stress  relationship  for  the  piezoelectric  gage. 

The  results  of  the  X-cut  LiNbO^  shots  show  that  there  is  a  large 
sensitivity  to  shear  stress;  however,  there  is  also  a  large  sensitivity 
to  compressive  stress,  despite  the  shorting  configuration,  that  renders 
this  gage  unusable  under  complex  loading  conditions. 

The  remaining  five  experiments  were  conducted  on  163°  Y-cut 
I.iNbO^  crystals.  These  gages  were  fabricated  to  our  specifications  and 
were  5  mm  thick.  The  use  of  thicker  gages  allows  the  use  of  thicker 
PMMA  specimens,  thereby  permitting  a  larger  temporal  separation  between 
the  compressive  and  shear  wave. 

The  results  of  shots  3  and  4  demonstrated  a  very  small  output 
from  the  compressive  wave;  however,  these  results  were  difficult  to 
interpret  because  of  large  impact  tilts.  The  impact  tilt  gives  rise 
to  the  following  complexities:  the  large  ratio  of  the  wave  velocity 
in  LiNbO^  to  the  projectile  velocity  results  in  a  large  inclination  of 
the  wavefront  and  an  additional  shear  wave  can  be  generated.  These 
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FIGURE  3  OUTPUT  FROM  X-CUT  LiNb03  GAGES  PLOTTED  AS  A  FUNCTION  OF 
NORMALIZED  TIME 

(a)  Compression  only.  Shot  *1,  (b)  Compression  and  Shear,  Shot  -2,  (c)  Difference 
of  Shot  *1  and  #2  to  show  signal  due  to  shear.  The  time  has  oeen  normalized 
with  respect  to  transit  time  through  the  gage. 


effects  and  the  lack  of  a  perfect  specific  direction  for  the  163° 

Y-cut  LiNbO^  give  rise  to  a  complex  wave  structure.  Rather  than 
analyze  this  wave-structure  (using  the  analysis  in  Section  II. B), 
we  chose  to  eliminate  the  problem  of  large  tilts. 

A  foam  flyer  plate  was  used  in  shots  5  and  6  to  permit  a  higher 
projectile  velocity,  thereby  reducing  the  closure  time  at  the  impact 
surface.  Particle  velocity  gages  were  vapour-deposited  on  the  Lucite 
specimen  plates  to  determine  the  input  stresses  in  the  Lucite.  Shot 
5  resulted  in  a  1.75-kbar  compressive  pulse,  and  a  negligible  signal 
output  from  the  gage  was  observed.  In  shot  6,  a  combined  compressive 
and  shear  experiment,  the  shear  particle  velocity  at  the  impact  surface 
was  measured.  Althcugh  the  signal  from  the  compressive  stress  was 
again  negligible,  the  shot  was  of  little  value  because  the  foam  imparted 
negligible  shear  motion  to  the  Lucite  specimen. 

In  shot  7  we  used  a  Lucite  flyer,  a  higher  projectile  velocity, 

and  a  larger  angle  (20°)  to  cause  larger  compressive  and  shear  stresses 

in  the  PMMA  and  hence  in  the  LiNbO^  gage.  An  assumed  linear  elastic 

response  in  compression  and  shear  for  the  PMMA  results  in  a  compressive 

stress  of  6  kbar  and  a  shear  stress  of  1  kbar,  respectively,  in  the 

* 

LiNbO^  gage.  However,  our  work  on  PMMA  shows  this  assumption  to  be 
wrong,  and  therefore  the  amplitude  of  the  shear  wave  in  PMMA  is  expected 
to  be  considerably  smaller  than  that  calculated  from  the  following 
linear  elastic  relation 


where 


T  =  p  CS  uS  (16) 

o  p  p 

g 

Cp  =  shear  wave  velocity 

5 

u^^  =  shear  particle  velocity  in  PMMA  at  the  impact 

surface. 


* 


We  are  currently  determining  the  shear  stress-strain  relation  in 
PMMA  under  a  separate  contract  with  the  U.S.  Army  Research  Office. 
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Figure  4  shows  the  voltage-time  results  from  shot  7  as  recorded  on  a 
dual-beam  oscilloscope.  The  upper  trace  is  the  LiNbO^  gage  output 
and  the  lower  trace  is  the  shear  particle  velocity  profile  from  the 
particle  velocity  gage  at  the  impact  surface.  The  time  settings  are 
0.2  ys/cm  for  both  traces;  the  voltage  settings  are  1  V/cm  for  the 
upper  trace  and  0.05  V/cm  for  the  lower  trace. 

The  particle  velocity  profile  at  the  impact  surface  (lower  trace) 
shows  that  initial  particle  velocity  imparted  to  the  specimen  is  90% 
of  the  value  calculated  by  assuming  a  "no  slip"  condition  at  impact. 

Only  the  first  0.8  ys  of  the  signal  is  useful  because  of  edge  rare¬ 
factions  that  influence  the  gage  leads.  The  decay  in  the  particle 
velocity  0.46  ys  after  impact  is  caused  by  the  shear  behavior  of  PMMA. 
This  decay  and  other  observations  suggest  that  PMMA  does  not  behave  in 
a  linear  elastic  manner,  and  the  LiNbO^  gage  output  from  shear  is 
expected  to  be  much  smaller  than  the  value  expected  for  a  1-kbar 
shear  stress. 

The  upper  trace  shows  a  small  signal  output  from  the  LiNbO^ 

gage  at  the  time  of  the  compressive  wave  arrival  at  the  gage.  The 

compressive  wave  velocity  in  PMMA  calculated  from  Figure  3  is  in 

excellent  agreement  with  independent  measurements  in  PMMA  made  at 
2  7 

SRI  and  elsewhere.  The  amplitude  of  the  signal  from  the  compressive 
wave  is  of  the  order  expected  from  the  small  piezoelectric  coefficient 
along  the  163°  Y-cut  direction.  In  future  work,  we  will  made  precise 
comparisons  between  the  theoretically  computed  values  and  the  experiment 
results. 

The  larger  signal  in  the  upper  trace  coincides  with  the  shear 

k 

wave  arrival  at  the  LiNbO^  gage.  This  result  thus  shows  that  the  163° 
Y-cut  LiNbO^  can  indeed  be  used  to  measure  the  shear  stress  at  the 
specimen-gage  interface.  As  explained  previously,  a  quantitative 
comparison  of  the  predicted  and  measured  gage  signal  is  not  possible 
at  present  because  of  lack  of  knowledge  of  the  dynamic  shear  strain 

k 

Shear  wave  velocity  in  PMMA  has  been  determined  from  other  independent 
measurements . 


MP-7306-4 


FIGURE  4  RESULTS  FROM  THE  163°  Y-CUT  LiNb03  IMPACT 
EXPERIMENT 

The  top  trace  represents  output  from  the  LiNb03 
gage  (1  V/div).  The  bottom  trace  represents  transverse 
particle  velocity  imparted  to  the  PMMA  (0.05  V/div). 
Time  scale  is  common  to  both  traces  (0.2  /isec/div). 
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behavior  of  PMMA  and  the  epoxy  bond.  The  decrease  in  signal  0.3  ys 
after  the  arrival  of  the  shear  wave  coincides  with  the  compressive 
wave  arriving  at  the  back  of  the  gage  and  represents  the  limit  of 
useful  recording  by  the  gage  in  our  experiment. 

The  results  from  the  163°  Y-cut  LiNbO^  experiments  may  be  summa¬ 
rized  as  follows:  these  gages  demonstrate  a  very  low  response  to 
compressive  stresses  and  a  large  response  to  shear  stresses  under 
dynamic  loading  conditions.  The  signal  from  the  shear  wave  is  nearly 
an  order  of  magnitude  higher  than  the  compressive  signal.  However, 
the  amplitude  of  the  compressive  wave  is  more  than  an  order  of  magnitude 
higher  than  the  shear  wave  amplitude.  These  results  confirm  our 
theoretical  calculations.  A  quantitative  calibration  of  the  gages 
requires  additional  experiments,  as  outlined  in  the  next  section. 


IV  DISCUSSIONS 


This  section  presents  a  discussion  of  the  theory  and  the  experi¬ 
mental  results,  as  well  as  specific  questions  that  need  to  be  answered 
for  the  development  of  a  gage  that  can  be  used  routinely  in  laboratory 
and  field  experiments. 

A.  Discussion  of  the  Experimental  Results 

The  results  of  the  two  experiments  on  X-cut  LiNbO^  clearly 
demonstrated  that  shear  sensitivity  of  the  crystal  alone  is  not  a 
sufficient  requirement  for  developing  a  gage  for  use  under  complex 
loading  situations.  Hence,  most  of  the  shear  transducers  commonly 
used  in  ultrasonics  cannot  be  used  for  our  purposes. 

A  fairly  general  analysis  of  alpha  quartz  and  LiNbO^  provided 

only  one  useful  gage  type:  163°  Y-cut  LiNb03,  which  satisfies  almost 

all  the  criteria  presented  in  Section  II. D.  The  shear-to-compressive 
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sensitivity  ratio  has  a  value  of  250.  Although  our  experiments  did 
not  provide  a  complete  quantitative  check,  they  did  demonstrate  the 
large  shear-to-compressive  sensitivity  ratio.  A  quantitative  comparison 
of  the  observed  signal  with  the  theoretical  calculations  can  be  made 
using  the  relation  between  current  and  stress  given  in  Appendix  A. 

In  future  work  we  will  carry  out  these  quantitative  comparisons. 
Furthermore,  in  future  analysis  of  the  experimental  results,  we  need 
to  consider  the  effect  of  the  small  deviation  of  the  163°  direction 
from  a  specific  direction. 

B.  Recommendations  for  Further  Gage  Development 

We  have  successfully  demonstrated  the  feafibility  of  developing 
a  piezoelectric  shear  stress  gage.  However,  the  routine  use  of  a  gage 
for  laboratory  and  field  applications  requires  further  assessment  of 
gage  response  accuracy,  quantitative  calibration,  and  proper  packaging 
(for  field  use).  These  considerations  are  briefly  discussed  below. 


The  primary  requirement  for  gage  development  was  the  determination 
of  the  appropriate  orientation.  The  orientation  calculations  were 
based  on  the  values  of  the  piezoelectric  constants.  For  LiNbO„, 
differences  in  the  constants  have  been  reported  in  the  literature. 

The  probable  cause  of  this  variation  is  the  material  source.  By 
performing  a  few  experiments  along  a  specific  direction  (165°  Y-cut), 
we  can  establish  the  accuracy  of  the  constants  without  introducing  the 
complexity  caused  by  the  deviation  from  specific  direction.  The  results 
of  these  experiments  will  provide  accuracy  bounds  for  gage  calibrations. 

In  using  the  piezoelectric  tensor  (Eq.  14  in  Section  II. D),  we 
need  to  examine  two  other  calibration-related  questions:  (a)  What  is 
the  effect  of  ^  0  on  the  gage  response?  This  arises  because  e^' 
is  not  zero  in  Equation  (14).  (b)  What  is  the  effect  of  shear  component 

being  non-zero  when  we  are  trying  to  measure  the  response  due  to 
0)?  The  second  question  is  important  in  field  conditions  when 
an  arbitrary  shear  strain  may  be  present  on  the  gage  face.  The  presence 
of  e'  can  give  rise  to  P'  ^  0. 

The  first  question  can  be  answered  by  conducting  precise  calibra¬ 
tion  experiments  that  would  also  establish  rhe  accuracy  of  the  gage 
response.  The  second  question  can  be  answered  by  performing  experi¬ 
ments  on  gages  that  are  rotated  90°  with  respect  to  gages  used  in  our 
work,  that  is,  e'  ^  0  and  e'  =  0  (Note,  e'  is  unchanged). 


Packaging  requirements  also  need  to  be  addressed  for  field  usage 
of  the  gage.  Our  examination  of  field  work  shows  that  measurements 
of  shear  stresses  at  soil-structure  interfaces  appear  to  be  most 
appropriate  with  the  piezoelectric  shear  stress  gage.  For  packaging, 
we  must  consider  optimal  placement  of  gages,  survivability,  and  tech¬ 


niques  for  attaching  leads  to  the  electrodes.  For  field  usage,  the 
gage  constants  are  best  expressed  in  terms  of  the  d...  ’s  and  care  must 
be  taken  to  ensure  the  proper  stress  boundary  conditions  (lateral 


stresses  should  be  zero) . 


Finally,  in  future  work,  we  should  also  examine  other  candidate 
materials.  Although  quartz  and  Lithium  Niobate  have  been  most 
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extensively  studied  for  large  amplitude  stress  measurements,  other 
materials  may  better  satisfy  the  criteria  outlined  in  Section  II. D 
for  a  shear  stress  gage. 
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Appendix  A 


OPERATION  MODES  FOR  PIEZOELECTRIC  GAGES" 


Piezoelectric  gages  may  be  used  in  either  of  two  modes:  current 
(or  short  circuit)  mode  and  charge  (or  open  circuit)  mode.  The  current 
mode  is  used  for  stress  measurements  for  time  durations  less  than  the 
transit  time  through  the  gage  (~  ys  time  scales) .  The  charge  mode 
is  used  for  stress  measurements  for  time  durations  (~  ms  time  scales) 
much  greater  than  transit  time  through  the  gage.  This  appendix  presents 
the  electrical  equations  for  each  of  these  modes  to  show  the  similarity 
in  the  basic  piezoelectric  concepts  for  gage  design.  The  assumptions 
used  in  deriving  these  relations  are  similar  to  those  used  by  Graham 
et  al."* 


Assuming  that  the  electric  fields  are  linear  and  one-dimensional, 
we  have 


D(h,t)  =  £(h,t)  *  E(h, t)  +  P(h,t)  (A.l) 


where  D  = 
e  = 

E  = 

P  = 
h  = 

If  we  assume 


Displacement  field 
Dielectric  permittivity' 

Electric  field 
Polarization  field 

Spatial  coordinate  ch rough  the  thickness 
constant  permittivity,  we  can  write 


o  o  o 


of  the  material. 


dh  (A. 2) 


Parts  of  the  material  presented  here  are  taken  from  a  review  by 
Y.  M.  Gupta,  "Response  and  Use  of  a  Quartz  Gage,"  Shock  Dynamics 
Laboratory  Internal  Report  73-03,  Washington  State  University  (1973). 

To  avoid  confusion  with  the  symbol  used  for  strains,  we  denote 
permittivity  with  a  bar  below  the  symbol. 


Jr r. 3£  oL^uUS-itCT 
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where  £  =  gage  thickness.  The  assumption  of  zero  conductivity  in  the 
gage  gives 


V  •  D  =  3D/9h  =  0 

Putting  this  expression  into  Eq.  (A. 2),  we  have 

I"  £  a 


(A.  3) 


Current  flow,  i,  between  the  gage  faces  may  be  expressed  using 
Maxwell's  equations 


*  -  -  -h  /  (?  • s  ] 


dv  =  A 


dD 

dt 


(A.  A) 


where  q  =  true  charge  density 

A  =  area  over  which  the  measurements  are  made. 

Eq.  (A. 4)  is  now  analyzed  with  respect  to  the  current  and  charge  mode 
usage. 

Current  Mode: 

The  short  circuit  condition  between  the  gage  faces  gives 

£ 


/ 


E  •  dh  =  0 


(A. 5) 


Combining  Eqs.  (A. 3),  (A. 4),  and  (A. 5),  we  have 

£ 


l  = 


_A_  _d_ 
£  dt 


/ 


P  •  dh 


(A. 6) 


If  we  assume  that  P  =  fa,  where  0  is  the  stress  component  of 

* 

interest  and  f  is  the  appropriate  piezoelectric  constant,  we  have 


This  scalar  relationship  is  only  possible  if  the  piezoelectric 
tensor  has  the  appropriate  form.  Details  of  relating  polarization 
fields  to  stresses  is  the  subject  of  Section  III.C. 
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i 


(A. 7) 


Because  the  gage  itself  is  assumed  to  be  linear  elastic,  we  can  write 


0(h,t)  =  o(h  -  U  ♦  t)  (A. 8) 

where  U  is  wave  velocity  in  the  gage.  Putting  Eq.  (A. 8)  in  (A. 7) 
and  recognizing  that  h  and  t  are  independent  variables,  we  obtain 


(o(£)  -  0(0)} 


(A. 9) 


Equation  (A. 9)  is  the  fundamental  equation  for  this  mode  and  relates 

the  current  output  to  stress  difference  between  the  two  faces  of  the 
gage.  For  times  less  than  the  transit  through  the  gage, 
o(£)  =  0  in  Eq.  (A. 9). 


Charge  Mode: 

The  open  circuit  condition  between  the  gage  faces  is  expressed 
by  setting  i  =  0  in  Eq.  (A. 4).  This  gives 


_d_ 

dt 


=  0 


(A. 10) 


Substituting  the  following  definitions  of  voltage  and  permittivity 


V 


J  E  dh  and  _£  = 


(A. 11) 


where  c  =  capacitance,  in  Eq.  (A. 10),  we  have 


V  = 


(A. 12) 
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Assuming  P  =  fa  as  before,  we  can  write 


V  = 


(A. 13) 


gage 


The  term  in  the  brackets  is  the  average  stress  o 
thickness. 


av 


through  the 


V 


f  AO 

av 

c 


(A. 14) 


Equation  (A. 14)  is  the  fundamental  equation  for  this  mode  and 
relates  the  voltage  output  to  the  average  stress  in  the  piezoelectric 
gage. 

Equations  (A. 9)  and  (A. 14)  show  that  both  modes  of  operation  use 
the  same  piezoelectric  relation  P  =  fa.  Hence,  the  piezoelectric 
developments  from  our  impact  work  (using  the  current  mode)  are  directly 
applicable  to  the  design  of  field  gages  that  would  use  the  charge  mode. 


Appendix  B 


TRANSFORMATION  OF  PIEZOELECTRIC  CONSTANTS 

This  appendix  describes  a  simple  and  convenient  method  to  determine 
the  crystallographic  orientations  best  suited  for  piezoelectric  gage 
applications.  Although  considerable  information  exists  on  piezoelectric 
constants  along  particular  directions,  the  method  presented  is  more 
general  and  convenient  for  the  applications  of  interest. 

The  third-rank  piezoelectric  tensors  defined  in  Section  XII  can 
be  transformed  as  follows 


=  a.  a.  a,  e 
ijk  iP  31  kr  Pqr 

(B.l) 

ijk  ai£  ajm  akn  ^lmn. 

(B.2) 

where  a  is  the  transformation  matrix  given  by  Eq.  (2)  in  Section  II. 
mn 

In  defining  the  coordinate  rotations  for  our  work,  we  have  chosen  the 
Institute  of  Radio  Engineers  (IRE)  conventions,  instead  of  the  usual 
solid  angles.  We  have  numerically  solved  equations  (B.l)  and  (B.2) 
for  coordinate  rotations  between  0°  and  180°  (in  increments  of  1°) 
about  the  X,  Y,  or  Z  axes.  The  numerical  program  has  the  following 
inputs. 

(1)  Reading  in  the  axis  of  rotation,  that  is,  1  or  2  or  3 
(Corresponding  to  X  or  Y  or  Z) . 

'2)  Reading  in  the  crystal  type  and  the  piezoelectric  matrix 
(numerical  values)  in  the  crystallographic  system. 

The  above  input  is  used  to  determine  the  rotated  constants  as 
f  olloxi/s 

(1)  Change  the  piezoelectric  constants  from  the  matrix  to  tensor 
notation. 

(2)  Construct  the  transformation  matrices  for  the  rotation  angle. 

(3)  Transform  the  piezoelectric  tensor. 
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(4)  Convert  the  transformed  tensor  back  to  the  matrix  form. 

(5)  Print  out  the  angle  of  rotation  and  the  transformed  piezo¬ 
electric  matrix. 

(6)  Repeat  the  above  steps  for  all  desired  angles. 

The  differences  in  converting  e  and  d  values  from  matrices  to 
tensors  and  vice-versa  are  included. 

As  an  example  of  this  procedure,  we  show  the  piezoelectric  stress 

matrix  (e..'s)  for  LiNbOn  rotated  about  the  X-axis  in  Table  B.l.  The 
ij  3 

values  shown  are  for  rotation  angles  between  0U  to  19°  and  152°  to  173°. 
The  non-zero  strains  in  the  impact  experiments  are  and  £^.  From 
the  matrices  in  Table  B.l  one  can  easily  determine  the  usefulness  of 
a  particular  orientation.  Note  the  marked  change  in  shear-to-compression 
sensitivity  ratio  (e^/e^)  '*'n  &oin8  ^rom  162°  to  164°. 


onn  own  own  of u  <*  owin  ow«  own  own  own  owo 

I  I  •»  l  »  l  l  I  i 
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